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ABSTRACT 

We develop a novel abundance matching method to construct a mock catalog of luminous red 
galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS), using catalogs of halos and subhalos 
C| in iV-body simulations for a A-dominated, cold dark matter model. Motivated by observations 

Q_f suggesting that LRGs are passively-evolving, massive early-type galaxies with a typical age 

^ 5 Gyr, we assume that simulated halos at z = 2 (z2-halo) are progenitors for LRG-host 
subhalos observed today, we label the most tightly bound particles in each progenitor z2- 
halo as LRG "stars". We then identify the subhalos containing these stars to z = 0.3 (SDSS 
redshift) in descending order of the masses of z2-halos until the comoving number density of 
the matched subhalos becomes comparable to the measured number density of SDSS LRGs, 
^lrg = 10~ 4 h 3 Mpc -3 . Our only free parameter is the number density of halos identified 
at z = 2 and this parameter is fixed to match the observed number density at z = 0.3. By 
tracing subsequent merging and assembly histories of each progenitor z2-halo, we can directly 
compute, from TV-body simulations, the distributions of central and satellite LRGs and their 
internal motions in each host halo at z = 0.3. While the SDSS LRGs are galaxies selected 
by the magnitude and color cuts from the SDSS images and are not necessarily a stellar- 
mass-selected sample, our mock catalog reproduces a host of SDSS measurements: the halo 
occupation distribution for central and satellite LRGs, the projected auto-correlation function 
of LRGs, the cross-correlation of LRGs with shapes of background galaxies (LRG-galaxy 
weak lensing), and the nonlinear redshift-space distortion effect, the Finger-of-God effect, 
in the angle-averaged, redshift-space power spectrum. The mock catalog generated based on 
our method can be used for removing or calibrating systematic errors in the cosmological 
interpretation of LRG clustering measurements as well as for understanding the nature of 
LRGs such as their formation and assembly histories. 

Key words: cosmology: theory - galaxy clustering - galaxy formation - cosmology: large- 
scale structure of the Universe 



1 INTRODUCTION 

Galaxy redshift surveys are one of the primary tools for studying 
the large-scale structure in the Universe (Davis & Huchra 1982; de 
Lapparent et al. 1986; Kirshner et al. 1987; York et al. 2000; Pea- 
cock et al. 2001). Over the coming decade, astronomers will have 
even larger surveys including BOSS 1 (Dawson et al. 2012), Wig- 
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gleZ 2 (Blake et al. 201 1), Vipers 3 , FMOS 4 , HETDEX 5 , BigBOSS e 
(Schlegel et al. 2009), Subaru Prime Focus Spectrograph (PFS 7 ; El- 
lis et al. 2012), Euclid 8 , and WFIRST 9 . The upcoming generation 
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of galaxy redshift surveys is aimed at understanding cosmic accel- 
eration as well as measuring the composition of the universe via 
measurements of both the geometry and the dynamics of structure 
formation (Wang et al. 1999; Eisenstein et al. 1999; Tegmark et al. 
2004; Cole et al. 2005). 

On large scales, galaxies trace the underlying distribution of 
dark matter, thus measurements of their correlations are an impor- 
tant probe of cosmology. Because of their relatively high spatial 
densities and their intrinsic luminosities, the luminous red galaxies 
(LRGs) are one of the most useful tracers (Eisenstein et al. 2001; 
Wake et al. 2006). Measurements of the clustering properties of 
LRGs have been used to measure the baryon acoustic oscillation 
(BAO) scale (Eisenstein et al. 2005; Percival et al. 2007; Ander- 
son et al. 2012) as well as to constrain cosmological parameters 
(Tegmark et al. 2004; Cole et al. 2005; Reid et al. 2010; Saito et al. 
2011). 

Our lack of a detailed understanding of the relationship be- 
tween galaxies and their host halos are the dominant systematic un- 
certainty in the analysis of large-scale clustering data. The halo oc- 
cupation distribution (HOD) approach or the halo model approach 
has been a useful, albeit empirical, approach to relating galaxies to 
dark matter (see e.g., Peacock & Smith 2000; Seljak 2000; Scoc- 
cimarro et al. 2001, for the pioneer works). In these approaches, 
the distribution of halos is first modeled for a given cosmologi- 
cal model, e.g. by using A-body simulations, and then galaxies 
of interest are populated into the dark matter halos. The previous 
works have shown that, by adjusting the model parameters, the 
HOD based model well reproduces the auto-correlation functions 
of LRGs measured from the Sloan Digital Sky Survey 10 (SDSS) 
(Zehavi et al. 2005; Zheng et al. 2007; Wake et al. 2008; Reid & 
Spergel 2009; White et al. 201 1). Based on these studies, it has been 
shown that LRGs reside in massive halos with a typical mass of a 
few times 10 13 h -1 Mq. However, the HOD method employs sev- 
eral simplified assumptions. For instance, the distribution of galax- 
ies is assumed to follow that of dark matter in their host halo and the 
model assumes a simple functional form for the HOD to constrain. 

In this paper, we pursue an alternative approach, abundance 
matching. The abundance matching method directly connects tar- 
get galaxies to simulated subhalos assuming a tight and physically- 
motivated relation between their properties, e.g., galaxy luminos- 
ity and subhalo circular velocity, without employing any free pa- 
rameters (e.g., Kravtsov et al. 2004; Conroy et al. 2006; Trujillo- 
Gomez et al. 2011; Reddick et al. 2012; Masaki et al. 2012; Nuza 
et al. 2012). However, it is not still clear whether the method can 
simultaneously reproduce different clustering measurements such 
as the auto-correlation function and the galaxy-galaxy weak lens- 
ing (Neistein & Khochfar 2012). Most authors only use the auto- 
correlation function to test their abundance matching model. 

However, rather than matching galaxies with halos at low red- 
shift, we take an alternative approach to abundance matching. Mo- 
tivated by observations suggesting that LRGs are passive, massive 
early-type galaxies, which are believed to have formed at z > 1 
(Masjedi et al. 2008; Carson & Nichol 2010; Tojeiro et al. 2012), 
we assume that the progenitor halos for LRG-host subhalos are 
formed at z = 2. (The conclusions of the method are not sensitive 
to the exact redshift choice). We identify the most massive halos 
at this redshift, follow the innermost particles of each progenitor 
halo, and then identify the subhalos at z = 0.3 containing these 
particles. We adjust the number of halos identified as LRG progen- 
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itors at z = 2 to match the observed number density of the SDSS 
LRGs, ft L RG ^ 10~ 4 h? Mpc" 3 (also see Conroy et al. 2008; Seo 
et al. 2008, for a similar idea based method). With this method, we 
can directly trace, from the simulations, how each progenitor halo 
at 2 = 2 experiences merger(s), is destroyed or survives at lower 
redshift as well as which progenitor halos become central or satel- 
lite galaxies in each host halo at z = 0.3. Thus, our method allows 
us to include assembly/merging histories of the LRG-progenitor 
halos. Our method is solely based on a mass-selected sample of 
progenitor halos at z = 2 and does not have any free parameter be- 
cause the mass threshold is fixed by matching to the number den- 
sity of SDSS LRGs. We compare statistical quantities computed 
from our mock catalog with the SDSS measurements: the HOD, 
the projected auto-correlation function of LRGs, the LRG-galaxy 
weak lensing and the redshift-space power spectrum of LRGs. Even 
though our method is rather simple, we show that our mock catalog 
remarkably well reproduces the different measurements simultane- 
ously. 

The structure of this paper is as follows. In Section 2, we de- 
scribe our method to generate a mock catalog of LRGs by using 
A-body simulations for a A-dominated cold dark matter (ACDM) 
model as well as the catalogs of halos and subhalos at z — 2 and 
z = 0.3. In Section 3, we show the model predictions on the rela- 
tion between LRGs and dark matter, and compare with the SDSS 
measurements. Section 4 is devoted to discussion and conclusion. 



2 METHODS 

2.1 Cosmological A-body simulations 

We perform two realizations of cosmological A-body simulations 
using the publicly-available Gadget-2 code(Springel et al. 2001b; 
Springel 2005). For each run, we assume a flat ACDM cosmology 
with Q. m = 0.272, Q b = 0.0441, fi A = 0.728, H = 100ft = 
70.2 kms -1 Mpc" 1 , a 8 = 0.807 and n 3 = 0.961 using the 
same parameters and notation as in the the WMAP 7-yr analysis 
(Komatsu et al. 2011). Our large scale simulation, which we la- 
bel at L1000, employs 1024 3 dark matter particles in a box of 
1 /i _1 Gpc on a side. This large box is most useful for statistical 
quantities such as correlation functions, so we will mainly use the 
LI 000 simulation to have enough statistics. We test convergence 
with a higher resolution simulation that employs 1024 3 particles in 
a box of 300 /i _1 Mpc on a side. This run is referred as L300. The 
mass resolution for the simulations (mass of an A-body particle) is 
7x 10 10 /i _1 M® or 1.9 x 10 9 h" x M Q for L1000 or L300, respec- 
tively. The initial conditions for both simulations are generated us- 
ing the second-order Lagrangian perturbation theory (Crocce et al. 
2006; Nishimichi et al. 2009) and an initial matter power spectrum 
at z = 65, computed from the CAMB code (Lewis et al. 2000). We 
set the gravitational softening parameter to be 30 and 8 /i _1 kpc for 
LI 000 and L300 runs, respectively. 

We use the friends-of-friends (FoF) group finder (e.g., Davis 
et al. 1985) with a linking length of b = 0.2 in units of the mean 
interparticle spacing to create a catalog of halos from the simula- 
tion output and use the SubFind algorithm (Springel et al. 2001a) 
to identify subhalos within each halo. In this paper, we use halos 
and subhalos that contain more than 20 particles. Each particle in 
a halo region is assigned to either a smooth component of the par- 
ent halo or to a subhalo, where the smooth component contains the 
majority of A-body particles in the halo region. Hereafter we re- 
fer the smooth component as the central subhalo and the subhalo(s) 
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as satellite subhalo(s), respectively. For each subhalo, we estimate 
its mass by counting the bounded particles, which we refer as the 
subhalo mass (M su b). We store the data (positions and velocities) 
of particles in halos and subhalos at different redshifts. To estimate 
the virial mass (M V j r ) for each parent halo, we apply the spher- 
ical overdensity method to the FoF halo, where the boundary of 
the spherical region is determined by the interior virial overdenisty, 
Avir, relative to the mean mass density (Bryan & Norman 1998). 
The overdensity A v i r — 268 at 2 = 0.3 for the assumed cosmolog- 
ical model. The virial radius is estimated from the estimated mass 
as .Rvir = (3M v ir/47rp m0 Avir) 1 ' /3 , where p m o is the comoving 
matter density. 

2.2 Mock catalog of LRGs: connecting halos at z = 2 to 
central and satellite subhalos at z — 0.3 

LRGs are very useful tracers of large-scale structure as they can 
reach a higher redshift, thereby enabling to cover a larger volume 
with a spectroscopic survey of LRGs for cosmology (Eisenstein 
et al. 2001, 2005). LRGs are passively-evolving, early-type mas- 
sive galaxies, and their typical ages are estimated as ~ 5 Gyrs 
(Kauffmann 1996; Wake et al. 2006; Masjedi et al. 2008; Carson & 
Nichol 2010). This implies that LRGs, at least a majority of the stel- 
lar components, were formed at z £ 1 (Masjedi et al. 2008). Moti- 
vated by this fact, we here propose a simplest abundance-matching 
method for connecting LRGs to dark matter distribution in large- 
scale structure as follows. 

Our method rests on an assumption that progenitor halos for 
LRG-host subhalos today are formed at z = 2, which is closer to 
the peak redshift of cosmic star formation rate (Hopkins & Beacom 
2006). Our choice of z = 2 is just the first attempt, and a formation 
redshift can be further explored so as to have a better agreement 
with the SDSS measurements (see Section 4 for further discus- 
sion). First, we select halos from the simulation output at z = 2 
as candidates of the progenitor halos (hereafter sometimes referred 
as 22-halo). In doing this, we select the 22-halos in descending or- 
der of their masses (from more massive to less massive) until the 
comoving number density becomes close to that of SDSS LRGs at 
2 = 0.3, which we set to tjlrg = 10 -4 /i 3 Mpc -3 . More exactly 
speaking, we need to identify more halos having the number den- 
sity of ~ 1.3 x 10~ 4 /i 3 Mpc~ 3 at least, because about 30% of 
22-halos, preferentially in a massive halo region at z = 0.3, ex- 
perience mergers from z — 2 to z = 0.3 for the assumed ACDM 
model (see below for details). Second, we trace the 30% innermost 
particles of each 22-halo particles to lower-redshift simulation out- 
puts until 2 = 0.3, where the innermost particles are defined by 
particles within a spherical boundary around the mass peak of a 
22-halo as can be found from Figure 1. Third, we make a match- 
ing of the traced innermost particles of each 22-halo to central and 
satellite subhalos at z = 0.3 (hereafter 20.3-subhalo). If more than 
50% of the innermost particles are contained in a 20.3-subhalo, 
we define the subhalo as the subhalo hosting LRG at z = 0.3. 
We repeat this procedure in descending order of the FoF masses 
of 22-halos until the comoving number density of the matched 
20.3-subhalos (LRG-host subhalos) is equal to the target value, 
mlrg = 10~ 4 /i 3 Mpc -3 . The minimum mass of LRG-progenitor 
halos at z — 2 is about 6 x 10 12 h~ x MQ (which contains about 
90 A-body particles for the LI 000 run). Thus we should stress that 
our method does not have any free parameter to adjust; the mass 
threshold for the LRG-progenitor halos is fixed by the target num- 
ber density. In addition, central and satellite subhalos are populated 
with galaxies under a single criterion: if a subhalo at z = 0.3 is a 



descendant of the 22-halo, the subhalo is included in the matched 
sample. On the other hand, the standard abundance matching uses 
different quantities for central and satellite subhalos, e.g., the cir- 
cular velocities at the redshift of target galaxies or at the accretion 
epoch, a epoch for a subhalo accreting onto the main host halo. In 
this way, we construct a mock catalog of LRGs at z = 0.3, and 
keep all the necessary information; the mass of LRG-progenitor 
22-halo, the mass of LRG halo at z — 0.3, the mass and location 
of each LRG-host subhalo and the internal motion inside the parent 
halo. With this method, we can trace the assembly history of each 
LRG-host halo as well as which LRGs become central or satellite 
galaxies in their parent halos at z = 0.3. We refer an LRG, which 
resides in a central or satellite subhalo, as a central or satellite LRG, 
respectively, in each parent halo. 

Some of the LRG-host halos at z = 0.3, especially mas- 
sive halos, contain multiple LRG-host subhalos in our mock cat- 
alog (see the example in the lower panel of Figure 1). We often 
call such systems "multiple-LRG systems" in the following dis- 
cussion (also see Reid & Spergel 2009; Hikage et al. 2012a). We 
refer the LRG-host halos, which host only a single LRG inside, 
as "single-LRG systems". The average halo masses for the single- 
and multiple-LRG systems are found from the LI 000 run to be 
M vir = 4.8 x 10 13 and 1.5 x 10 14 /i _1 M 0> respectively. The frac- 
tion of the multiple-LRG systems among all the LRG-host halos is 
about 8% in the L1000 run. Assuming that most stars in each LRG 
are formed until 2 = 2 and the total stellar mass scales with host 
halo mass, i.e. mass of 22-halo, we will often refer the most massive 
22-halo in a multiple-LRG system as the brightest LRG (BLRG), 
while we refer the smallest 22-halo as the faintest LRG (FLRG). 
Note that we also refer an LRG in a single-LRG system as BLRG. 
A BLRG does not necessarily become a central galaxy in the par- 
ent halo at 2 = 0.3. Similarly, a central LRG in a multiple-LRG 
system is not necessarily the most massive 22-halo, i.e. BLRG, al- 
though the central subhalo is the most massive subhalo in the parent 
halo by definition. 

Table 1 summarizes properties of the LRG-host halos in our 
mock catalog, computed using the L1000 run. In the table, we give 
the average virial mass and radius for all the LRG-host halos and 
for each sample of the single- and multiple-LRG halos. The table 
also gives the fraction of satellite LRGs or central LRGs in each 
sample of the host halos. These results can be compared with the 
recent results found in Hikage et al. (2012a) from the SDSS DR7 
LRG catalog, as given in the table caption. The mock catalog fairly 
well reproduces the SDSS results. In the following, we study these 
quantities in more detail. 

Figure 1 shows snapshots of the A-body particle distribution 
in the L1000 run outputs at different redshifts, for the regions where 
multiple- or single-LRG systems are formed at z = 0.3. The figure 
clearly shows how each LRG-progenitor halo is defined at z = 2, 
how the innermost particles are traced to lower redshift and how 
LRG-progenitor halos merge with each other and become to reside 
in central and satellite subhalos at the final redshift z = 0.3. Our 
method allows us to directly trace the merging and assembly his- 
tories of LRG-progenitor halos. Although the number density of 
LRG-host subhalos is set to the density of LRGs as we described 
above, the figure shows that more LRG-progenitor halos or subha- 
los survive (with less mergers) at higher redshift than at z = 0.3. 
Hence our method also has a capability to study what kinds of halos 
or subhalos at higher redshift are progenitors for the SDSS LRGs 
(see Section 4 for a further discussion). 

Figure 2 shows how each LRG-progenitor halo at z = 2 loses 
or acquires its mass due to merging and/or tidal stripping when it 
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Type of LRG-host halos 


Total number 


M vir [r^g] 


R vh . [ft^Mpc] 


/sat — LRG 


BLRG 
Veen 


„FLRG 
ycen 


All LRG-host halos 


91,090 


5.64 x 10 13 


0.80 


0.099 


0.96 




Single-LRG halos 


83,891 


4.81 x 10 13 


0.78 


0.022 


0.98 




Multiple-LRG halos 


7,199 


1.52 x 10 14 


1.1 


1.00 


0.74 


0.21 



Table 1. Summary of properties of LRG-host halos, computed from the mock LRG catalog in the L1000 simulation run (see text for details). Here we consider 
all LRG-host halos and the single- and multiple-LRG halos that host only one or multiple LRG(s) inside, respectively. M v i r and R v i r are the average virial 
mass and radius of the host halos. / sa t-LRG is the fraction of halos that have satellite LRG(s) among all LRG-host halos taken in each row. <7 B = n RG is the 
fraction of halos that host its BLRG as a central galaxy among all the halos, where BLRG is the brightest LRG, the most massive LRG-progenitor halo at 
2 = 2, compared to other LRG-subhalo(s) in the same host halo at z = 0.3. Note that, for the single-LRG hosts, we call the LRG as the BLRG. g^ c „ RG is 
the fraction of halos that host FLRG as a central LRG, where FLRG is the faintest LRG, the smallest LRG-progenitor halo, in each multiple-LRG halo. Note 
that the fraction of satellite LRGs among all the LRGs (not the LRG-host halos) in the mock catalog is 0.108. The results can be compared with the results in 
Hikage et al. (2012a) for the SDSS DR7 LRG catalog: the fraction of multiple-LRG systems among all the LRG-host halos is 0.045, the average halo masses 
for the single- and multiple-LRG systems are M v ir — 0.36 and 1.36 X 10 14 h~ 1 MQ, respectively, and the fractions of central BLRGs and FLRGs in the 
multiple-LRG systems g^ C n RG ~ 0.54 and iJccn RG — °- 32 - respectively. 




Figure 1. Evolution of dark matter (./V-body) particle distribution around the region of subhalos hosting mock LRGs at 2 = 0.3, taken from our LI 000 
simulation run. The upper-row panels are for the region around the host-halo of the brightest LRG among single-LRG systems (the host halo mass M v j r = 
8.42 X 10 14 h~ 1 Mq), systems which host only one LRG inside in the z = 0.3 output, while the lower-row panels are the most massive host-halo among 
systems hosting one central and three satellite LRGs (Mvir = 1.44 X 10 15 /i _1 Mg). The dot symbols in each panel are member particles in the halo regions 
at 2 = 2 or the subhalo region(s) at lower redshifts. The red-color particles are the 30% innermost particles of each halo at z = 2 and selected based on our 
abundance matching method between the progenitor halos and the LRG-host subhalos at z = 0.3 to reproduce the observed number density of SDSS LRGs 
(see text for details). Then we trace where the red-color particles are distributed in each subhalo region at lower redshift. By matching the red-color particles 
to central and satellite subhalos in each host halo of z = 0.3 output, we can define locations of each LRG in a host halo at z = 0.3; if a subhalo at z = 0.3 
contains more than 50% of the red-color particles of a progenitor halo, we define it as an LRG-host subhalo. The upper-row panels show the case that 1 1 
progenitor halos of LRGs are formed at z = 2, and then are merged at lower redshift, forming one central LRG in the host halo at 2 = 0.3. The lower-row 
panels show that 24 progenitor halos at 2 = 2 form one central LRG and three satellite LRGs in the host-halo at 2 = 0.3. The blue circles in the panel of 

1/3 

2 = 0.3 shows the positions of mock LRGs. The size of each circle is proportional to M , , where M su t, is the subhalo mass. 
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Figure 2. Compai'ison between masses of the LRG-progenitor halos at z = 2 and the LRG-host subhalos at z = 0.3, computed from the L1000 ran, where 
each progenitor halo and subhalo are matched based on our method (see Figure 1). The left and right panels show the results for all the LRG-host halos and the 
multiple-LRG systems, respectively. The black and red points are for central and satellite LRGs, respectively. Note that the central LRG-subhalo is a smooth 
component of the parent halo at z = 0.3. The line in each panel shows the case that the progenitor halo does not either acquire or lose its mass until z = 0.3: 
Msub( z = 0.3) = Mp p(z = 2). The figure shows that satellite LRGs preferentially reside in less massive progenitor halos at z = 2, some subhalos for 
satellite LRGs lose their masses due to tidal stripping when accreting into more massive halos, and subhalos for central LRGs acquire their masses due to 
merger. The upper- and right-side panels in each plot are the projected distributions for central and satellite LRGs along the y- or rr-axis direction, respectively. 



becomes an LRG-host subhalo at z = 0.3, computed using the 
catalogs of halos and subhalos in the z — 2 and z — 0.3 outputs of 
LI 000 run. Note that the halo mass shown in the x-axis, Mf f, is 
the FoF mass, the sum of FoF particles in each halo region. First, 
the figure shows that we need to select the LRG-progenitor halos at 
z = 2 down to a mass scale of about 6 x 10 1 h~ Mq. It is shown 
that some subhalos for satellite LRGs lose their masses due to tidal 
stripping as implied in Figure 1, while subhalos for central LRGs 
acquire their masses due to merger. Comparing the left and right 
panels manifests that multiple-LRG systems tend to reside in more 
massive LRG-progenitor halos at z = 2 and become more massive 
LRG-host halos at z = 0.3, and that the mass difference between 
subhalos for central and satellite LRGs is larger in multiple-LRG 
systems, implying a larger difference between their luminosities 
(see Hikage et al. 2012a, for a similar discussion). 

Thus our method is solely based on the masses of LRG- 
progenitor halos at z = 2 (see Figure 2) and the connection with 
central and satellite subhalos in the parent halos at z = 0.3. On 
the other hand, LRGs in the SDSS catalog are selected based on 
the magnitude and color cuts from the SDSS imaging data (pri- 
marily gri), and are not necessarily a stellar-mass-selected sample, 
although their stellar masses are believed to have a tight relation 
with the host halo masses. Nevertheless, we will show below that 
the mock catalog surprisingly well reproduces the SDSS measure- 
ments. 

Since LRGs in our mock catalog reside in relatively massive 
halos at z = 2, with masses MfoF £ 6 x 10 12 h~ 1 M Q (Figure 2), 
as well as in massive parent halos at z — 0.3, our method does 
not necessarily require a high spatial resolution for simulations. A 
simulation with 1 /j Gpc on a side length seems sufficient, which 
allows an accurate estimation of statistical quantities of LRGs. This 
is not the case if one wants to work on the abundance matching 



method for less massive galaxies or more general types of galaxies 
(e.g., Trujillo-Gomez et al. 201 1 ; Reddick et al. 2012; Masaki et al. 
2012). 



3 RESULTS: COMPARISON WITH THE SDSS LRG 
MEASUREMENTS 

3.1 Halo occupation distribution and properties of satellite 
LRGs 

First, we study the halo occupation distribution (HOD) for LRGs in 
Figure 3, where the HOD gives the average number of LRGs that 
the halos at z — 0.3 host as a function of host-halo mass. Here we 
consider the HODs for central and satellite LRGs which reside in 
central and satellite subhalos in the LRG-host halos, respectively. 
Again we should emphasize that our method does not assume any 
functional forms for the HODs, unlike done in the standard HOD 
method, and rather allows us to directly compute the HODs from 
the simulations. The figure shows that, even if LRG-progenitor ha- 
los are selected from halos at z = 2 by a sharp mass threshold, 
our mock catalog naturally predicts that the central HOD has a 
smoother shape around a minimum halo mass, as a result of their 
merging and assembly histories from z — 2 to z — 0.3. To be more 
precise, the central HOD is smaller than unity (HOD < 1) for host 
halos with M vil £ 10 14 /i _1 M s , meaning that only some fraction 
of the halos can host a central LRG. The more massive halos can 
host multiple LRGs inside. Conversely, the L1000 run shows that 
the fraction of massive halos, which do not host any LRG, is 1.3% 
for halos with masses M vir ^ 1 x 10 14 /i -1 M , while all halos 
with Mvir > 2 x 10 14 fe _1 M have at least one LRG inside. 

In Figure 3, we compare the HODs of our mock catalog 
with the result in Reid & Spergel (2009, hereafter RS09), where 
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Figure 3. The halo occupation distribution (HOD) for LRGs as a function 
of parent halo mass, measured from our mock catalog. Our mock catalog 
has an assignment of each LRG to central or satellite subhalos in a parent 
halo at z = 0.3 (see Figure 1), thereby allowing us to directly compute the 
HODs for central (solid curve) and satellite (dashed) LRGs from the mock 
catalogs. The black and blue curves are the results from the L1000 and 
L300 runs, respectively, where the L300 run is a higher resolution ran with a 
small box size, 300 /r'Mpc (see text for details). The red curves show the 
SDSS measurements, taken from Reid & Spergel (2009, RS09). RS09 fixed 
the function form of central HOD, and then constrained the satellite HOD 
from the SDSS LRG catalog using the Counts-in-Cylinders technique. The 
hatched region is the range allowed by varying each model parameter of 
the satellite HOD within its lcr confidence range. The mock catalog well 
reproduces the SDSS measurements, including the shape of central HOD 
around the cutoff mass scale as well as the slope and amplitude of satellite 
HOD, without employing any free parameter to adjust after the abundance 
matching. 
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Figure 4. The fraction of halos hosting satellite LRG(s) inside as a function 
of halo mass, computed by using all the LRG-host halos at z = 0.3 in the 
L1000 andL300 runs. 



the HOD was constrained by using the Counts-in-Cylinders (CiC) 
method for identifying multiple LRG systems from the SDSS DR7 
LRG catalog with the aid of halo catalogs in TV-body simulations. 
Although the cosmological models are slightly different in between 
our study and RS09 and the assumed redshift is slightly different 
(z = 0.2 in their study), we employed the same best-fit parameters 
to compute the LRG HOD in RS09 for this figure. 

Due to limited constraints from the SDSS LRG catalog, espe- 
cially for low-mass host-halos, RS09 assumed the fixed form for 
the central HOD: 



(iVcen(M)) 



1 



1 + erf 



log M - log Afmin 
Clog M 



(1) 



with M m i n = 5.7 x 10 h Mo and a losM = 0.7, in order 
to obtain meaningful constraints on the satellite HOD. The central 
HOD for low-mass host-halos is difficult to constrain, because low- 
mass host-halos of LRGs are observationally difficult to identify. 
Our mock catalog naturally predicts a smoother shape for halos 
with M < 10 14 /i _1 M Q for the central HOD, which fairly well 
agrees with the assumption in RS09. Therefore, we do not think 
that the difference for the central HODs is significant, and needs to 
be further carefully studied. 

On the other hand, the satellite HOD in RS09 is almost per- 
fectly recovered by our mock catalog, where RS09 assumed the 
functional form for the satellite HOD to be given by (JV S at (AT)) = 
(iVcen(M)) [(M - M cu t)/Mi] a and then constrained the param- 
eters (M cut , Mi, a) from the SDSS LRG catalog. The hatched re- 
gion is the range at each host-halo mass bin that is allowed by vary- 
ing the model parameters within the lcr confidence regions. Our 
results confirm that parent halos of ~ 10 15 cj _1 A/q have up to 
several LRGs inside, as first pointed out in RS09. The L300 result, 
the simulation result with higher spatial resolution, gives similar 
results to the L1000 results, showing that the numerical resolution 
is not an issue for the HOD results. Even though SDSS LRGs are 
selected by the magnitude and color cut, not by their host subhalo 
masses, our method seems to capture the origin of SDSS LRGs; 
mass-selected halos at z ~ 2 are main progenitors of LRGs, and 
their subsequent assembly and merging histories determine where 
LRGs are distributed within the host halos at lower redshift. 

One motivation of this paper is to understand the physics of the 
nonlinear redshift-space distortion, i.e. the Finger-of-God (FoG) ef- 
fect, in the redshift-space power spectrum of LRGs. The FoG effect 
is caused by internal motion of satellite LRG(s) in LRG-host halos 
(Hikage et al. 2012b,a). In the following, we study several quanti- 
ties relevant for the FoG effect; the fraction of satellite LRGs, the 
radial profile of satellite LRGs inside the parent halos and the inter- 
nal velocities of satellite LRGs (see Hikage et al. 2012b, for details 
of the theoretical modeling). 

Figure 4 shows the fraction of LRG-host halos at 2 = 0.3 
hosting satellite LRG(s) inside, as a function of halo mass. Note 
that we include the single-LRG systems, which host one LRG as a 
satellite galaxy, in addition to the multiple LRG systems. The er- 
ror bars around the solid curve are Poisson errors, estimated using 
the number of halos in each mass bin. The figure shows that more 
massive halos have a higher probability to host satellite LRG(s). 
About 20% of parent halos with M vil ~ 10 14 h~ 1 M G host satel- 
lite LRG(s). 

We naively expect that BLRG, the most massive LRG- 
progenitor halo at z = 2 among LRG-progenitor halo(s) in the 
same LRG-host halo at z = 0.3, becomes a central galaxy in the 
host halo. The solid curves in Figure 5 show the fraction of BLRGs 
for becoming a satellite galaxy in LRG-host halos at z = 0.3 as 
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Figure 5. The solid curves show the fraction of the parent halos hosting the 
brightest LRG (BLRG) as a satellite galaxy, among all the LRG-host halos. 
Here the BLRG is the most massive LRG-progenitor halo at z = 2 among 
all the progenitor halos which become to reside in the same LRG-host halo 
at z = 0.3. The dashed curves are the similar fraction of LRG-host halos 
with satellite BLRG, but computed using only the multiple-LRG systems. 
The error bars are computed from the number of halos in each mass bin 
assuming Poisson statistics. 
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Figure 6. The average radial profile of satellite LRG host subhalos, ob- 
tained by stacking the positions of satellite LRGs in all the LRG-host halos 
with satellite LRG(s), as a function of radius relative to the virial radius of 
each parent halo. The mean mass of the LRG-halos used in this calculation 
isAf vir ~ 1. 31 orl.24xl0 14 /i- 1 M Q for theLlOOO or L300 runs, while 
the mean virial radius is R vlI ~ 1.07 or 1.06 h _1 Mpc, respectively. For 
comparison, the upper dotted curve shows the profile of dark matter aver- 
aged for the same host halos with an arbitrary amplitude. The error bars 
at each radial bin are estimated by first dividing LRG-host halos into 27 
subsamples (27 subvolumes) and then computing the variances at each ra- 
dial bin. The typical off-center radius for satellite LRGs appears to be about 
400 /i -1 kpc. 
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Figure 7. The average radial velocity of satellite LRGs, (v B,r\ with re- 
spect to the halo center in each LRG-host halo, computed by using all the 
LRG-host halos with satellite LRG(s) as in the previous figure. The neg- 
ative (v s tr \ means a coherent infall towards the halo center. The upper 
curves show the average radial velocity dispersion around the coherent in- 
fall, cr of j ,.. For the comparison, the dotted curve shows the average ve- 
locity dispersion expected from virial theorem, <j v - lr = v/ GM vlr /2R v i r . 

The combination of (v f[ ir \ and o" D ff r implies that satellite LRG(s) sink 
towards the halo center due to dynamical friction, and then have an os- 
cillating motion around the halo center with the velocity dispersion of 
~ 500 km s" 1 . 



a function of halo mass, computed using all the LRG-host halos. 
For parent halos with M vir > 10 14 /i _1 M , there is up to 10% 
probability for its BLRG to be a satellite galaxy. 

The dashed curves are the similar fraction, but computed us- 
ing only the multiple LRG systems. This sample is intended to 
compare with the recent result in Hikage et al. (2012a) (see be- 
low). In this case, the fraction of satellite BLRGs is higher for host 
halos with smaller masses, with larger error bars. This can be ex- 
plained as follows. Recall that most of low-mass host-halos with 
masses 10 14 /i~ 1 A/q are single-LRG systems as can be found 
from Figure 3. That is, only a small number of the halos become 
multiple-LRG systems, causing larger Poisson error bars at each 
mass bin. We have found from the simulation outputs that such 
low-mass halos of multiple LRG systems (mostly the systems with 
2 LRGs) tend to display a bimodal mass distribution due to on- 
going or past major merger, where the BLRG and other (mostly 
central) LRG tend to have the small mass difference. As a result, 
such low-mass multiple-LRG systems have a higher chance to host 
the BLRG as a satellite LRG. On the other hand, the fraction of 
the halos hosting satellite BLRG converge to the solid curve with 
increasing the host-halo mass, because most of such massive ha- 
los are multiple-LRG systems. For multiple LRG systems with 
M vir ~ 10 14 /i _1 M , about 30% of BLRGs are satellite galax- 
ies. 

Recently, Hikage et al. (2012a) studied the multiple-LRG sys- 
tems defined from the SDSS DR7 catalog by applying the CiC 
technique as well as the FoF group finder method to the dis- 
tribution of LRGs in redshift space. Then they used the differ- 
ent correlation measurements, the redshift-space power spectrum, 
the LRG-galaxy lensing and the cross-correlation of LRGs with 
photometric galaxies, to study properties of satellite LRGs. They 
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found that the multiple-LRG systems has a typical halo mass of 
Misoi, ~ 1.6 x 10 14 /i _1 M (with a roughly 10% statistical 
error), and that about 40% of BLRGs in the multiple-LRG sys- 
tems appear to be satellite galaxies 11 . Using the method in Ap- 
pendix C in Hu & Kravtsov (2003), we find M vir ~ 0.85Mi 80 b 
assuming an NFW profile with concentration parameter c = 4. 
The weak lensing for the SDSS multiple-LRG systems implies 
M v i r = 1.4 x 10 14 7i _1 M from Misof, = 1-6 x 10 14 /i _1 M . 
Our mock catalog shows a fairly good agreement with the SDSS re- 
sults, for the average halo mass and the fraction of satellite BLRGs 
(also see Table 1). 

In Figure 6, we study the average radial profile of satellite 
LRGs. In this calculation, we employ only the parent halos con- 
taining satellite LRG(s), and estimate the radial profile by stacking 
the radial distribution of satellite LRG(s) in units of the radius rel- 
ative to the virial radius of each halo. We use the mass peak of the 
smooth component as the halo center. The average profile p a g is 
normalized as 

j dr 47rr' 2 p off (r') = 1, with r' = r off /i? vir , (2) 

where r ff is the distance from the density maximum of the smooth 
component. The average mass of the host halos is M v i r — 1.31 x 
10 14 or 1.24 x 10 14 7i _1 M for the L1000 or L300 run, respec- 
tively, while the average virial radius _R v ir — 1-07 or 1.06 /i -1 Mpc 
in the comoving unit. Compared to the dark matter profile, the ra- 
dial profile of satellite LRGs clearly displays a flattened profile. The 
typical off-center radius, where the profile starts to be flattened, is 
found to be about 400 /i kpc as 7? v j r ~ 1 h~ 1 Mpc, which is in 
a good agreement with the result in Hikage et al. (2012a). The ra- 
dial profile also shows a decline at the smaller radii. Thus our result 
clearly shows a violation of the assumption often used in a standard 
HOD method that the radial profile of member galaxies follows the 
dark matter profile. However, note that the L300 run shows no satel- 
lite LRG at small radii r g/R v i T ^ 0.1, except for the innermost 
bin. The satellite LRGs at the small radii are found from the most 
massive host-halos in LI 000 run, which do not exist in the smaller 
box simulation, L300. Although the mock catalogs show a sharp 
rise at the innermost bin r a s/R vlI ~ 0.06 (~ 60 /i _1 kpc), which 
may indicate merging LRGs to the central subhalo, the scatters are 
large even for the L300 run, so the result is not significant. Never- 
theless, it is worth mentioning that the satellite LRG distribution in 
our mock catalog seems to show a similar profile to the profile of 
most massive subhalos in cluster-scale halos in Gao et al. (2012) 
(see Figures 15 and 16 for the profile). These features in the radial 
profile of massive subhalos may be as a result of dynamical fric- 
tion, tidal stripping and merging to the central galaxy. However, this 
needs a further careful study to derive a more robust conclusion, by 
using high-resolution simulations as well as a larger number of the 
realizations. 

Figure 7 shows the average radial profile of internal motions 
of satellite LRGs in the parent halos, where the bulk motion of 
each parent halo (the average velocity of iV-body particles belong- 
ing to the smooth component of the halo) is subtracted from the 
velocity of each LRG-host subhalo. We considered only the host 
halos with satellite LRG(s) as in Figure 6. The curves, labelled as 
(v s,r), are the average radial velocities for all the satellite LRGs 

11 Note that Hikage et al. (2012a) used the halo mass definition of Mi8o(>, 
which is different from the virial mass M v ; r and defined by the enclosed 
mass inside which the mean density is 180 times the mean background mass 
density. 
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Figure 8. Projected auto-correlation function of LRGs, w p (R), as a func- 
tion of the projected distance R. The solid and dashed curves show the 
results from our mock catalogs in the LI 000 and L300 runs, respectively. 
The error bars are estimated using the measurements from 8 subdivided vol- 
umes of each simulation volume. The square and diamond symbols are the 
correlation functions measured from the SDSS catalog of LRGs at z ~ 0.3, 
taken from Zehavi et al. (2005) and Masjedi et al. (2006), respectively. 



with respect to the halo center. The average velocity is negative, 
reflecting the coherent infall motion towards the halo center, and 
the infall velocity is larger with increasing radius up to the virial 
radius. The average radial velocity becomes zero on average at the 
halo center. These support that the LRG-host subhalo approaches 
to the halo center due to dynamical friction. On the other hand, the 
curves, labelled as a a B,r, are the average velocity dispersions of 
satellite LRGs. The velocity dispersion has greater amplitudes with 
decreasing the radius, reaching to a g >r ~ 500 km s _1 . For com- 
parison, the horizontal dotted line shows the average virial velocity 
dispersion, <7 v ir = \J GM v \ r /2R v \ r = 521 kms -1 , for the ha- 
los hosting satellite LRG(s) in the LI 000 run. The combination of 
{VoS,r) and a g ir implies that satellite LRGs gradually approach 
to the halo center due to dynamical friction and have an oscillating 
motion around the halo center. Again the amplitude of the velocity 
dispersion, a Q s, r — 500 km s _1 , is in nice agreement with the re- 
cent measurement in Hikage et al. (2012a), where they found the 
velocity dispersion of about 500 km s _1 for satellite LRGs in the 
multiple-LRG systems. 



3.2 Projected correlation function 

Next we study the projected auto-correlation function of LRGs, 
w p (R), defined as 

/•7T max 

w p (R) =2 d-n £ 99 (r = a/^T~R2). (3) 

Jo 

where R is the projected separation between two LRGs in the pairs 
used for the correlation measurement in units of the comoving 
scale, 7T is the separation parallel to the line-of-sight and £ 99 (r) is 
the three-dimensional correlation function. Following Zehavi et al. 
(2005), 7r ma x is set to be 80 h~ Mpc. The projected correlation 
function is not affected by the redshift-space distortion effect due 
to peculiar velocities of LRGs. 
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Figure 9. The average surface mass density profile around LRGs, which 
is an observable of the LRG-galaxy weak lensing. The solid and dashed 
curves are the results of our mock catalog, obtained by stacking TV-body 
particles around all the LRG-host subhalos in the LI 000 and L300 runs, 
respectively. The error bars are estimated using the measurements from 27 
subsamples of LRG-host subhalos. The data with error bars show the SDSS 
measurements in Mandelbaum et al. (2012). 

In Figure 8, we compare the projected correlation function 
measured from our LRG mock catalog with the SDSS measure- 
ments (Zehavi et al. 2005; Masjedi et al. 2006). In the SDSS mea- 
surements, Zehavi et al. (2005) used an LRG sample in the magni- 
tude range of —23.2 < M g < —21.2 and with the mean redshift 
(z) ~ 0.3. Masjedi et al. (2006) used the same sample to extend the 
measurement down to very small scale, below R — 500 /i _1 kpc, 
by taking into account various observational effects such as the 
fiber collision. Note that the number density of the SDSS LRGs 
is 9.7 x 10 -5 /i 3 Mpc -3 , which is slightly smaller than the number 
density of the mock LRGs (10 -4 /i 3 Mpc -3 ), and the cosmologi- 
cal model employed in the measurement is slightly different from 
the model we assumed for our simulations. The figure shows that 
our mock catalog remarkably well reproduces the projected cor- 
relation function of LRGs over a wide range of separation radii, 
including the 1- and 2-halo terms that arise from correlations be- 
tween LRGs within the same host halo and in different host halos, 
respectively. Comparing the results for the L1000 and L300 runs 
reveals that the correlation function for LI 000 has a smaller ampli- 
tude at R < 0.7 h~ 1 Mpc than that for L300. Thus the L1000 run 
implies a systematic error due to the lack of numerical resolution at 
the small scales. The L300 result shows a better agreement with the 
SDSS measurement in Masjedi et al. (2006). The small-scale clus- 
tering arises mostly from multiple-LRG systems, different LRGs 
in the same halo, so that numerical resolution seems important to 
resolve these small subhalos. 



3.3 LRG-galaxy weak lensing 

Correlating the positions of LRGs with shapes of background 
galaxies, the so-called LRG-galaxy weak lensing, is a powerful 
means of probing the average dark matter distribution around the 
LRGs (Mandelbaum et al. 2006, 2012). The LRG-galaxy lensing 
measures the radial profile of differential surface mass density de- 
fined as 



AT,(R) = £(< R) - E(7?). (4) 

The profile E(_R) is the average surface mass density around the 
LRGs defined as 

E(-R) = p m0 J dn[l + £ gm (r = y/^ + R 2 )], (5) 

where p m o is the mean background mass density today, and £ 9m (r) 
is the three-dimensional cross-correlation between LRGs and the 
surrounding matter. In Eq. (4), £(< R) is the surface mass density 
averaged within a circular aperture of a radius R. Our use of the 
mean mass density today (p m o) is due to our use of the comoving 
units. 

Figure 9 shows that the average mass profile measured for all 
LRGs in the mock catalog is in good agreement with the SDSS 
measurement in Mandelbaum et al. (2012). Note that, to obtain the 
average mass profile from our mock catalogs, we stacked all N- 
body particles around all the LRG-host subhalo in the simulations, 
including the particles outside dark matter halos. The lensing signal 
at the radii smaller than about 1 h~ Mpc arises from the mass 
distribution within the same halo, while the signal at the larger scale 
arises from the mass distribution surrounding the host halos - the 
so-called 1- and 2-halo terms (e.g. see Oguri & Takada 201 1). The 
mock catalog well reproduces the measurements of different scales. 

Hikage et al. (2012a) also used the SDSS LRG catalog to 
study the weak lensing for the multiple-LRG systems. When mak- 
ing the lensing measurements, they used three different proxies 
for the halo center of each multiple-LRG system, the BLRG, 
FLRG and the arithmetic mean position of member LRGs (here- 
after "Mean"). By comparing the lensing signals for the different 
centers, they constrained the average off-center profiles for satellite 
BLRGs and FLRGs in the parent halos, finding about 400 /i _1 kpc 
for a typical offset radius. Figure 10 shows that the mock catalog 
predictions are in remarkably good agreement with the SDSS mea- 
surements for the different centers. Since these lensing signals are 
from the same parent halos of multiple-LRG systems, the differ- 
ences between the signals of different centers arise from the off- 
centering effects. As nicely shown in Hikage et al. (2012a), the 
lensing signals for the BLRG and FLRG centers can be well ex- 
plained by a mixture of the central and satellite LRGs in the sam- 
ple. The lensing signals for the FLRG center have smaller am- 
plitudes due to a larger fraction of satellite (off-centered) FLRGs 
than in the BLRG centers. On the other hand, the Mean center 
does not have any galaxy (subhalo) at its position, and therefore 
the Mean center always has an off-centering effect from the true 
center in each LRG system. This causes decreasing powers of the 
lensing signal at the smaller radii than the typical off-center radius. 
The lensing signals at some radii for the FLRG and Mean centers 
show some discrepancy from the mock catalog, but we do not think 
that the disagreement is significant. The average masses inferred 
from the SDSS measurement and the mock catalog (Table 1) for 
the multiple-LRG halos agree within about 30%; M v ir — 1-36 or 
1.52 x 10 14 h~ 1 M Q , respectively. 

As can be shown in Figures 8, 9 and 10, our mock catalog of 
LRGs well reproduces both the SDSS measurements for the auto- 
correlation function of LRGs and the LRG-galaxy weak lensing si- 
multaneously. As recently discussed in Neistein & Khochfar (2012) 
(also see Neistein et al. 201 1), the abundance matching method has 
a difficulty to reproduce these measurements with the same model, 
although they considered the spectroscopic sample of SDSS galax- 
ies, rather than focused on LRGs. Thus the agreements of our mock 
catalog show a capability of our method to predict different statis- 
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Figure 10. The average surface mass profiles for the multiple-LRG systems. The different panels show the results obtained by taking the different centers in 
each multiple-LRG halo; the brightest LRG (BLRG), the faintest LRG (FLRG) and the center-of-mass of different LRGs or the arithmetic mean positions of 
member LRGs (Mean) in the left, middle and right panels, respectively. The data with error bars show the SDSS measurements for the multiple-LRG systems 
in Hikage et al. (2012a). 



tical quantities of LRGs by self-consistently modeling, rather than 
assuming, the fractions of satellite LRGs among different halos and 
the radial distribution of satellite LRGs in the host halos (also see 
Masaki et al. 2012, for a recent development on the extended abun- 
dance matching method based on the similar motivation). 

3.4 Redshift-space power spectrum of LRGs 

Another observable we consider is the redshift-space power spec- 
trum of LRGs. The FoG effect due to internal motion of galaxies 
is one of systematic errors to complicate the cosmological interpre- 
tation of the measured power spectrum. The FoG effect involves 
complicated physics inherent in the evolution and assembly pro- 
cesses of galaxies, so is very difficult to accurately model from the 
first principles. The standard way to reduce the FoG contamination 
is to remove satellite galaxies from the region of each multiple- 
LRG system, and to keep only one galaxy (LRG in our case), ide- 
ally the central galaxy, because the central galaxy is supposed to 
be at rest with respect to the parent halo center and does not suf- 
fer from the FoG effect (see below for the related discussion). For 
example, Reid et al. (2010) developed a useful method for this pur- 
pose; first, reconstruct the distribution of halos from the measured 
distribution of LRGs by identifying multiple-LRG systems based 
on the CiC and FoF group finder method, and then keep only one 
LRG for each multiple-LRG system. However, the chosen LRG 
is not necessarily the central galaxy (more exactly, they used, as 
the halo center proxy, the arithmetic mean of member LRGs or 
the center-of-mass of different CiC groups without any luminosity 
or mass weighting), so there may generally remain a residual FoG 
contamination in the measured LRG power spectrum as pointed out 
in Hikage et al. (2012a). 

In the left panel of Figure 11, we study the FoG effect on 
the redshift-space power spectrum, caused by the off-centering ef- 
fect of LRGs in our mock catalog. Following the method in Reid 
et al. (2010) and Hikage et al. (2012a), we study the redshift-space 
power spectrum for halos that host LRG(s), instead of the power 
spectrum for LRGs. To compute the power spectrum of halos, we 
need to specify the halo center in each LRG-host halo. For single- 
LRG systems, we use the LRG position as the halo center proxy. 



For multiple-LRG systems, we employ different proxies of halo 
center for each system as done in Figure 10; BLRG, FLRG or 
the arithmetic mean (Mean), where the Mean center is computed 
in redshift space taking into account redshift space distortion due 
to peculiar velocities of LRG-host subhalos. The figure shows the 
angle- averaged redshift-space power spectra for the different cen- 
ters, relative to the power spectrum for the mass center of each 
LRG-host halo (the mass center of TV-body particles of the host 
halo). Note that, for the power spectrum measurement, we used the 
exactly same number of LRG-host halos, and the difference is only 
the positions of halo centers in the multiple-LRG systems. Hence, 
the difference between the different spectra should be from the 
off-centering effects in the multiple-LRG systems. The spectra for 
BLRG, FLRG and Mean centers all show smaller amplitudes with 
increasing wavenumber, as expected in the FoG effect. To be more 
precise, the power spectrum of FLRG center shows the strongest 
FoG effect, because a larger fraction of FLRGs are satellite galax- 
ies than BLRGs (see Table 1). These results can be compared with 
Figure 2 in Hikage et al. (2012a). The mock catalog is in a qualita- 
tively good agreement with the SDSS measurements: the spectra of 
BLRG and Mean centers are similar, and the spectrum for FLRG 
shows the stronger FoG suppression. 

Figure 1 1 also shows the spectrum when taking the potential 
minimum as the center of each halo, where the potential minimum 
is the mass density peak of smooth component, i.e. the central sub- 
halo position, in each LRG-host halo. This corresponds to the case 
that the power spectrum is measured by using the position of a cen- 
tral galaxy in each host halo. (BLRG is not necessarily a central 
galaxy as shown in Figure 5). The power spectrum for the poten- 
tial minimum has a smaller amplitude than that of the mass center 
of host halo, implying that the potential minimum has a random 
motion around the mass center in each halo. Comparing the spec- 
tra for the potential minimum and the BLRG center shows that the 
BLRG spectrum has a smaller amplitude than the spectra for the 
potential minimum or the mass center by a few % in the fractional 
amplitude up to k ~ 0.3 h Mpc -1 . The few %-level FoG contami- 
nation would be okay for a current accuracy of the power spectrum 
measurement, but will need to be carefully taken into account for 
a higher-precision measurement of upcoming redshift surveys. For 
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Figure 11. The angle-averaged redshift-space power spectra for the LRG-host halos at z = 0.3, computed from the L1000 run. The different curves show 
the fractional differences of the power spectra using different proxies of each LRG-host halo position in the power spectrum estimation, relative to the power 
spectrum using the mass center as the halo position. Left panel: The dotted, dashed and dot-dashed curves are the results when using different halo center 
proxies for each multiple-LRG system; the arithmetic mean position of the member LRGs in redshift space (Mean), BLRG or FLRG as in Figure 10. Note 
that we use the LRG position as the halo center for each single-LRG system. Thus the differences between the different spectra arise only from the different 
positions of multiple-LRG systems in redshift space, to be compared with Hikage et al. (2012a). The different power spectra show decreasing amplitudes 
with increasing wavenumber, which is caused by the nonlinear redshift-space distortion, the so-called Finger-of-God effect, due to the internal motions of the 
chosen halo centers in LRG-host halos. We also show the power spectrum measured using the potential minimum of each LRG-host halo, where the potential 
minimum is the mass density peak of the smooth component of the halo that is likely to host the central galaxy. For comparison, the three dots-dashed curve 
shows the effect caused by massive neutrinos assuming the total neutrino mass m^tot = 0.104 eV. Right panel: Similar to the left panel, but the power 
spectrum using all the LRGs is added (the three dots-dashed curve). The power spectrum includes contributions from multiple LRGs in the same halo. To 
have a fair comparison, we properly subtracted the shot noise due to the different number densities of the LRG-host halos or the LRGs, 1/«lrg- halos ° r 
l/«LRGs- 
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Figure 12. Similarly to the previous figure, but for the halved samples of LRG-host halos. Left panel: The LRG-halos are divided into two halved samples by 
the halo masses; one sub-sample is defined by halos which have masses smaller than the median mass ("small-half"), while the other is by halos with masses 
larger than the median ("massive-half"). The massive halo sub-sample shows a stronger FoG effect. Right panel: Similar plot, but using only the single-LRG 
halos. 



comparison, the three dots-dashed curve shows the effect on the 
power spectrum caused by massive neutrinos, where we assumed 
m^tot — 0.1 eV for the total mass of neutrinos, close to the lower 
limit on the neutrino mass for the inverted mass hierarchy. For the 
normal mass hierarchy, the lower limit on the total mass is about 



0.05 eV, and the amount of the suppression is about half of the re- 
sult of 0.1 eV in Figure 11. 

In the right panel of Figure 11, we also show the redshift- 
space power spectrum if using all the LRGs in the catalog. Note 
that we properly subtracted the shot noise contamination from the 
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measured power spectra by using the number densities of LRGs 
or LRG-host halos. In this case, the power spectrum ratio shows 
greater amplitudes with increasing wavenumber rather than the 
FoG suppression. That is, the LRG power spectrum shows a greater 
clustering power or greater bias than in the LRG-host halo spec- 
trum. Since the scales shown here are greater than a few tens 
Mpc and therefore much greater than a virial radius of most mas- 
sive host-halos, the 1-halo term arising from clustering between 
two LRGs in the same host-halo should not be significant at these 
scales. Hence, the greater amplitudes in the LRG power spectrum 
would be due to a more weight on more massive halos, which have 
larger biases, because satellite LRGs preferentially reside in more 
massive halos. Since the effect of different linear bias should cause 
only a scale-independent change in the power spectrum ratio, the 
change in the LRG power spectrum should be from a stronger non- 
linear bias of such massive halos, overcoming the stronger FoG 
suppression in such massive halos. In fact, a combination of the per- 
turbation theory of structure formation and halo bias model seems 
to reproduce such a non-trivial behavior in the power spectrum am- 
plitudes (Nishizawa et al., in prep.). The results in the figure imply 
that including satellite LRGs in the power spectrum analysis com- 
plicates the interpretation of the measured power spectrum. 

In Figure 12, we study how the residual FoG effect varies with 
masses of LRG-host halos. To study this, we divide the LRG ha- 
los into two sub-samples by masses of the LRG-halos smaller and 
larger than the median, and measured the fractional power spectra 
for each sub-sample relative to the halo sample. As expected, the 
FoG effect is larger for the sub-sample containing more massive ha- 
los, because of the higher fraction of satellite BLRGs as well as the 
larger velocity dispersion (larger halo mass). The right panel shows 
the similar results, but obtained only by using the single-LRG ha- 
los. First of all, the single LRG systems have a smaller FoG effect, 
because of the smaller fraction of satellite BLRGs (Figure 5 and Ta- 
ble 1) as well as the smaller velocity dispersions for the lower-mass 
host-halos. Among the single-LRG halos, more massive halos have 
relatively a larger FoG contamination, but only by a few percent at 
k ^ 0.35 h Mpc -1 in the amplitude. Thus, the use of single-LRG 
systems may allow a cleaner interpretation of the measured power 
spectrum, yielding a more robust, unbiased constraint on cosmo- 
logical parameters. 



4 DISCUSSION AND CONCLUSION 

In this paper, we developed a new abundance-matching based 
method to generate a mock catalog of the SDSS LRGs, using cat- 
alogs of halos and subhalos in iV-body simulations. A brief sum- 
mary of our method is as follows: (1) identify LRG-progenitor ha- 
los at z = 2 down to a certain mass threshold until the comoving 
number density of the halos become similar to that of the SDSS 
LRGs atz — 0.3 (2) trace the merging and assembly histories of the 
LRG "star particles", the 30% innermost particles of the halo mem- 
ber particles that lie closest to the center of each LRG-progenitor 
halo, and (3) at z = 0.3, identify the central and satellite subhalos 
hosting the LRG "star" particles. If a subhalo at z = 0.3 contains 
more than 50% of the innermost particle of the progenitor halo, we 
assigned the halo at z = 2 and the subhalo at z = 0.3 as the LRG- 
progenitor halo and the LRG-host subhalo, respectively. We should 
emphasize that our method does not employ any free parameter 
to adjust in order for the model to match the measurements, once 
the mass threshold of the LRG-progenitor halos is determined to 
match the number density of SDSS LRGs. Thus, by assuming that 



a majority of stellar components of LRG is formed at z = 2, we 
can trace the assembly and merging histories of LRGs over a range 
of redshift, z = [0.3, 2]; for example, we can directly trace which 
LRGs become central or satellite galaxies in the LRG-host halos 
at z = 0.3. The novel aspect of our approach is that the abun- 
dance matching of halos to a particular type of galaxies (LRGs in 
this paper) is done by connecting the halos and subhalos at differ- 
ent redshifts (z = 2 and z = 0.3 in our case), while the standard 
method is done for the same or similar redshift to the redshift of 
target galaxies. In addition, central and satellite subhalos are popu- 
lated with LRGs under a single criterion: if a subhalo at z = 0.3 is 
a descendant of the z2-halos, the subhalo is included. The standard 
abundance matching uses different quantities for central and satel- 
lite subhalos, e.g., the circular velocities at the galaxy redshift and 
at the accretion epoch, respectively. 

Using the mock catalog, we computed various statistical quan- 
tities: the halo occupation distribution, the projected correlation 
function of LRGs, the mean surface mass density profile around 
LRGs (which is an observable of the LRG-galaxy weak lensing), 
and the redshift-space power spectrum of LRGs. We showed that 
the mock catalog predictions are in a good agreement with the mea- 
surements from the SDSS LRG catalog (Figures 3, 8,9, 10 and 11). 
Thus our method seems to capture an essential feature of LRG for- 
mation in terms of a hierarchical structure formation scenario of 
ACDM model. 

In the SDSS sample, only 5% of the halos contain multiple 
LRGs. In our simulation, 8% of the halos contain multiple LRGs. 
This modest deviation may be due to one of our many simple as- 
sumptions. First, we assumed that LRG progenitors are formed at 
a single epoch, z = 2. This is too simplified assumption, and a 
formation of LRGs should have some duration. Including some 
duration of formation of LRG-progenitor halos may improve the 
model prediction. Second, although LRGs are observationally se- 
lected by the well-designed magnitude and color cuts and not nec- 
essarily a mass-selected sample in terms of the stellar masses, our 
definition of the LRG-progenitor halos at z = 2 is solely based 
on their masses. The agreements between our mock catalog and 
the SDSS measurements support that the matching based on the 
LRG-progenitor halo masses seems fairly reasonable to mimic a 
population of LRGs. However, the model may be further refined 
by combining masses of the progenitor halos with other indicators 
when matching to LRGs. For instance, using the maximum circu- 
lar velocity of each halo instead of its mass may improve the model 
accuracy. Third, we did not use satellite subhalos in the parent halo 
at z = 2 in making the abundance matching. We naively expect 
that subhalos at z = 2 merge into central subhalos from z = 2 to 
z = 0.3 due to dynamical friction, so we used the simplest method 
as the first attempt. However, including the subhalos at z = 2 for 
the abundance matching may improve an accuracy of the mock cat- 
alog. These would be worth exploring. 

Our mock catalog or more generally our abundance matching 
method offers several applications to measurements. First, Masjedi 
et al. (2008) showed that, by using the small-scale clustering signal 
and the pair counting statistics, LRGs are growing by about 1.7% 
per Gyr, on average from merger activity from z — 1 to z ~ 0.3. 
Our method allows us to directly trace how each LRG-progenitor 
halo acquires the mass from other LRG-progenitor halos by major 
or minor mergers from z — 2 to z = 0.3. Hence, we can compare 
the prediction of our mock catalog with the measurement for the 
mass growth rate of LRGs. By using the constraint, we may be able 
to further improve the mock catalog. 

Second, our method can predict how the distribution of LRG- 
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progenitor evolves in relative to dark matter distribution as a func- 
tion of redshift. Thus, our mock catalog can be used to predict vari- 
ous cross-correlations of LRG positions with other tracers of large- 
scale structure. As one such example, in this paper, we studied the 
LRG-galaxy weak lensing measured via cross-correlation of LRGs 
with shapes of background galaxies, and showed a remarkably 
good agreement between our model and the SDSS measurements. 
Another cross-correlation that has been studied in the literature is a 
cross-correlation of LRGs with a map of cosmic microwave back- 
ground (CMB) anisotropies, which probes the stacked Sunyaev- 
Zel'dovich (SZ) effect (Hand et al. 201 1; Sehgal et al. 2012) or the 
lensing effect on the CMB. Since every massive halos always host 
at least one LRG (100% of halos with M vil )2x 10 14 /i _1 M Q 
in our mock catalog), the cross-correlation is a powerful cross- 
check of the SZ signals. Our mock catalog can predict how the 
stacked SZ signals change for different catalogs of LRGs such as 
an inclusion of satellite LRGs and multiple-LRG systems, which 
may be able to resolve some tension between the observed LRG- 
CMB cross-correlation signal and the theoretical expectation (Seh- 
gal et al. 2012). 

Third is an application of our method to LRGs or massive 
red galaxies at higher redshift than z = 0.3. The SDSS-III BOSS 
survey is now carrying out an even more massive redshift survey 
of SDSS imaging galaxies. The magnitude and color cuts used 
for the BOSS survey are designed to efficiently select galaxies at 
0.4 < z < 0.7, and are different from the SDSS-I/II LRG selec- 
tion. The BOSS galaxies are called the "constant mass" (CMASS) 
galaxies. The majority of CMASS galaxies are early-type galax- 
ies, but are not exactly the same population as LRGs. In addition, 
the comoving number density of CMASS galaxies is higher than 
that of LRGs by a factor of 3. Hence, it would be interesting to 
apply the method developed in this paper to the CMASS galax- 
ies. Figure 1 shows an interesting indication of our mock catalog: 
more LRG-progenitor halos survive in the z — 0.5 output than at 
z — 0.3, because the halos do not have enough time to experience 
merging due to the shorter time duration from z — 2. Hence, our 
mock catalog naturally predicts a higher number density of LRG- 
progenitor halos at higher redshift than at z — 0.3, and may be able 
to match some of the BOSS galaxies without any fine tuning. Since 
the BOSS survey will provide us with a higher-precision cluster- 
ing measurement and therefore has the potential to achieve tighter 
cosmological constraints, it is critically important to use an accu- 
rate mock catalog of the CMASS galaxies in order to remove or 
calibrate various systematic errors inherent in an unknown relation 
between the CMASS galaxies and dark matter. We hope that our 
method is useful for this purpose and can be used to attain the full 
potential of the BOSS survey or more generally upcoming redshift 
surveys for precision cosmology. This is our future study and will 
be presented elsewhere. 
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